Perturbative QCD Calculations of Elliptic Flow and Shear Viscosity 
in Au+Au Collisions at ^snn = 200 GeV 
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The elliptic flow V2 and the ratio of the shear viscosity over the entropy density, rj/s, of gluon 
matter are calculated from the perturbative QCD (pQCD) based parton cascade Boltzmann ap- 
proach of multiparton scatterings. For Au+Au collisions at -^s — 200A GeV the gluon plasma 
generates large 172 values measured at the BNL Relativistic Heavy Ion Collider. Standard pQCD 
yields rj/s ~ 0.08 — 0.15 as small as the lower bound found from the anti-de Sitter/conformal field 
theory conjecture. 
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The large values of the elliptic flow V2 measured at 
the Relativistic Heavy Ion Collider (RHIC) [ij indicate 
that the matter created, the quark gluon plasma (QGP), 
behaves as a nearly perfect fluid exhibiting strong "ex- 
plosive" collective motion. As transport coefficients like 
the shear viscosity characterize the dynamics of the QGP 
, it is obvious that a less viscous medium will generate 
larger V2 values. Now, ideal hydrodynamical calculations 
have reproduced the observed V2 dependence on central- 
ity (within 30%), mass, and transverse momentum px 
(for pT < 1.5 GeV) 0]. Hence, the QGP created pos- 
sesses a small viscosity coefficient [1]. The smallness of 
the viscosity is of great interest as a result of the recent 
debate about speculative "realizations" of supersymmet- 
ric representations of Yang-Mills theories using the anti- 
de Sitter/conformal field theory (AdS/CFT) conjecture 
0. On the other hand, calculations done within ideal hy- 
drodynamics failed to reproduce the V2 at higher central- 
ity, at higher transverse momentum, and its dependence 
on rapidity Q. Although recent calculation within vis- 
cous hydrodynamics 0] can fit the centrality dependence 
of V2 , the pt dependence cannot be obtained consistently. 
One drawback in these calculations is the assumed ideal 
thermal initial conditions. In principle, the entire process 
of the nuclear collision is theoretically better studied by 
using a consistent microscopic transport theory. Trans- 
port simulations can provide the correct viscous hydrody- 
namics when the system is near equilibrium. Moreover, 
nonequilibrium processes, like equilibration and the on- 
set and the breakdown of viscous hydrodynamics, are 
involved. However, both hadronic Q and partonic 
cascade calculations with 2^2 processes failed to yield 
the large V2 values observed at RHIC. 

Recently we have developed a relativistic perturbative 
QCD (pQCD) based on-shell parton cascade Boltzmann 
approach of multiparton scatterings (BAMPS) [lol . 11 1 
to study, on a semiclassical level, the dynamics of gluon 
matter produced in Au+Au collisions at RHIC energies. 
Although coherent quantum effects like color instabili- 



ties 12| may play a role at the very initial stage where 
the matter is superdense, it is the pQCD interaction that 
drives the gluon matter to equilibrium within a time in- 
terval of 1 fm/c. The fast thermalization observed is due 
to gluon bremsstrahlung [ll[. This process is also re- 
sponsible for the small ratio of the shear viscosity to the 
entropy density, 77/3 , for a static gluon gas [l3| . For a cou- 
pling constant as = 0.6(0.3) we obtain 77/s ~ 0.076(0.13), 
which is as small as the lower bound ry/s = l/47r found 
from the AdS/CFT conjecture ^. Hence, pQCD inter- 
actions can explain why the QGP created at RHIC is 
strongly coupled and behaves like a nearly perfect fluid. 
There is no need to invoke exotic black hole physics in 
higher dimensions and supersymmetric Yang-Mills theo- 
ries using the AdS/CFT conjecture to understand RHIC 
data. 

In this Letter we employ BAMPS to calculate the ellip- 
tic flow V2 of gluons in Au+Au collisions at ^/s = 200A 
GeV. For a firm footing we compare our results with the 
experimental data, assuming parton-hadron duality. To 
see how viscous the gluon plasma behaves, the ratio of 
shear viscosity to entropy density is calculated. 

The initial gluon distributions are taken as given in 



Ref. ll|: an ensemble of minijets with transverse mo- 



menta greater than 1.4 GeV, produced via semihard 
nucleon-nucleon collisions. We use Glauber geometry 
with a Woods-Saxon profile and assume independent bi- 
nary nucleon-nucleon collisions. A formation time for 
initial minijets is also included [lol 11|. In contrast 



to the Glauber-type initial condition, the color-glass- 
condensate-type initial condition provides larger initial 
eccentricity and thus might yield larger V2 values 14 1. 



Including quarks in BAMPS is in progress. We do not 
expect large changes in the present results shown below, 
because the quark amount at RHIC is initially only 20%. 
As quarks are less interactive than gluons, the total V2 
will be slightly smaller than that obtained in a pure gluon 
matter. 

Gluon interactions included in BAMPS are 
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elastic scatterings with the pQCD cross section 
da^s-^sg/dql = Qna^Jiql + "^^))^ as weh as pQCD 
inspired bremsstrahlung gg ggg with the effective 
matrix element 113 



12gWi_ 



2 (qi + m2,)2 ki[(ki-q^)2+^ 
Q{k±Ag — cosh y) , 



(1) 



where g^ — Airag. q± and kj^ denote the perpen- 
dicular component of the momentum transfer and of 
the radiated gluon momentum in the center-of-mass 
frame of the collision, respectively, y is the momen- 
tum rapidity of the radiated gluon in the center-of- 
mass frame, and Ag is the gluon mean free path, which 
is calculated self-consistently [13] ■ The interactions 
of the massless gluons are screened by a Debye mass 
mjj = ndcasNcJ (Pp/{2tt)^ ■ f/p, where dc = 16 is 
the gluon degeneracy factor for Nc = 3. tod is cal- 
culated locally using the gluon density function / ob- 
tained from the BAMPS simulation. The suppression 
of the bremsstrahlung due to the Landau-Pomeranchuk- 
Migdal effect is taken into account within the Bethe- 
Heitler regime using the step function in Eq. ([1]). 

In the present pQCD simulations the interactions of 
the gluons are stopped when the local energy density 
drops below 1 GeV/fm'^. This value is assumed to be 
the critical value for the occurrence of hadronization, 
below which parton dynamics is not valid. Because 
hadronization and hadronic cascade are not yet included 
in BAMPS, a gluon, which ceases to interact, propagates 
freely and can be regarded as a free pion according to the 
parton-hadron duality. 

In Table H] the initial and final transverse energy per ra- 
pidity dET/dy{\y\ < 1.5) for 6 centrality bins are listed, 
which are obtained from BAMPS using two different cou- 
pling constants. 

The PHOBOS values are extracted from the px spectra 
for charged hadrons [3], i.e., dEx = J dpT \/p'^ + mgX 
dN/dpT with Too being the pion mass. To compare with 
them we have multiplied a factor of 2/3 to dEx/dy for 
gluons. The dEx/dy using as — 0.3(0.6) is slightly 
larger (smaller) than the PHOBOS data within ±8%. 



TABLE I: Centrality dependence of the transverse energy 
per rapidity at midrapidity. Dimension is GeV. 

Centrality PHOBOS Initial Final = 0.3 Final = 0.6 
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Figure [T] shows the rapidity dependence of the ini- 
tial and final dEx/dy (solid curves with open symbols) 
for the most central 5% collisions, compared with the 
BRAHMS data for the charged meson rapidity distribu- 
tions, dEr/dy = niTdN/dy \]Jj- The shape of dEx/dy 
obtained from BAMPS agrees well with the BRAHMS 
data within y < 2.5. 

The chosen initial condition of minijets, with po = lA 
GeV as a parameter, seems appropriate for describing 
the experimental data. However, calculations with the 
same initial condition and employing only elastic pQCD 
00 gg interactions show much less decrease of dEx/dy 
[lO| (see dashed curves with solid symbols in Fig. [1]). 

Figure [2] shows the elliptic flow V2 at midrapidity and 
its dependence on rapidity, as obtained from the BAMPS 
calculations with as = 0.3, and 0.6. These are compared 
with the PHOBOS data [ISj]. 

Except for the central centrality region the results with 
as — 0.6 agree perfectly with the experimental data, 
whereas the results with ag = 0.3 are roughly 20% 
smaller. The results with only gg —f gg interactions 
(dashed curves with solid symbols) are much smaller than 
the data, which demonstrates that bremsstrahlung and 
backreactions are the dominant processes in generating 
large elliptic flow. 

Figure [3] shows the generation of the elliptic flow at 
midrapidity for noncentral Au-t-Au collisions with impact 
parameter 5 = 8.6 fm corresponding to (A^part) = 111. 

V2 increases strongly during the early stage of thermal 
equilibration (time scale ~ 1 fm/c) where hydrodynamics 
is not valid but where the eccentricity is largest. The gen- 
eration rate of V2 during thermalization is only slightly 
smaller than that during the subsequent (viscous) hydro- 
dynamical expansion. The V2 generation ends at 3 fm/c, 
where the critical energy density is reached for all the 
gluons at midrapidity. 

Hence, the generation of the large elliptic flow observed 
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FIG. 1: (color online) . Rapidity dependence of the transverse 
energy per rapidity. Results from BAMPS are multiplied by a 
factor of 2/3 to compare with the BRAHMS data for charged 
mesons. 
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FIG. 2: (color online). Upper panel: Elliptic flow «2(|yj < 1) 
from BAMPS using Us — 0.3 and 0.6, compared with the 
RHIC data [3. Lower panels: Pseudorapidity dependence 
of elliptic flow for three centrality classes ranging from pe- 
ripheral to central (25% - 50%, 15% - 25%, and 0% - 15%) 
from top to bottom. 



at RHIC is well described by pure perturbative gluon in- 
teractions as incorporated in BAMPS. It is consistent 
with the decrease of the transverse energy due to the 
mechanical work done by the pressure during the ex- 
pansion. Inclusion of hadronization [19^ and subsequent 
hadronic cascade ^.20] will yield only moderate contribu- 
tions to the final elliptic flow values, as the densities are 
already rather low and the eccentricity becomes small. 

The ratio of the shear viscosity to the entropy density, 
ij/s, is extracted from the BAMPS calculations (includ- 
ing gg <-> ggg) to see how viscous the gluon matter is for 
generating the large V2 values that match the experimen- 
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FIG. 3: (color online). Generation of elliptic flow at midra- 
pidity for a noncentral collision with an impact parameter of 
b = 8.6 fm. 



tal data. We obtain the shear viscosity as [IS 
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(2) 
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where n is the gluon density, E is the gluon energy, and 
Vz — Pz/E is the velocity in the beam direction. X)-^*' 
denotes the total transport collision rate [11] describing 
momentum isotropization. R23 and R32 are interaction 
rates for gg — > ggg and its backreaction, respectively. 
The entropy density is extracted by assuming that the 
gluon matter is in local kinetic equilibrium, which gives 
s = An — n In A, where A = n/ricq defines the gluon fugac- 
ity. The true entropy density is expected to be (slightly) 
smaller, because overall kinetic equilibration cannot be 
complete in an expanding system. 

Figure 3] depicts the shear viscosity over entropy den- 
sity ratio r]/s extracted locally at the central region dur- 
ing the entire expansion before decoupling. 

The central region is the region with transverse ra- 
dius of 1.5 fm and within the space-time rapidity interval 
\rit\ < 0.2, where ?7t = 5 ln[(t + z)/{t — z)]. The curves in 
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FIG. 4: (color online) . Shear viscosity to entropy density ratio 
rj/s at the central region. The upper band shows the results 
with Oa = 0.3 and the lower band the results with Oa = 0.6. 
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Fig|4]are obtained from BAMPS simulations with impact 
parameter & = 4, 6.3, and 8.6 fm, for as = 0.3 and 0.6. 
Note that during thermal equilibration (time scale ~ 1 
fm/c) Eq. (21) for calculating the shear viscosity is not 
fully valid and the true entropy density will be smaller 
than that given by the formula used. 

The r]/s ratio obtained from the simulations with the 
same as and for various impact parameters are almost 
the same, although the gluon density with b — 8.6 fm is 
significantly smaller than that with 6 = 4 fm. Interaction 
rates and transport collision rates scale with the temper- 
ature T. The gluon fugacity is close to 1. Hence, rj/s 
depends practically only on as- This also explains that 
the rj/s values extracted in the outer transverse regions 
(not shown) are nearly the same as the ratio in the center. 
Figure H] shows rj/s ~ 0.15 for as = 0.3 and rj/s ~ 0.08 
for as = 0.6. These values agree well with those found 
for a static thermal gluon gas [13 1. 

Figure [5] shows the cross sections of gg — > gg and gg 
ggg processes, extracted in the central region. 



E 



=0.3, b=4 fm 
=0.3, b=6.3fm 
=0.3, b=8.6 fm" 
=0.6, b=4 fm 
=0.6, b=6.3fm- 
=0.6, b=8.6fm 





t (fm/c) 



t (fm/c) 



FIG. 5: (color online). Cross sections of gg — > gg and gg 
ggg processes in the central region. 

Indeed, the cross section becomes larger when the den- 
sity is smaller, so that (cr) depends practically only on 
as- Note that the cross sections are indeed small. This 
demonstrates that it is not necessary to use unphysically 
large cross sections or ultrashort hadronization times 
[§| to yield large V2 values. 

The pQCD based parton cascade BAMPS is used to 
calculate the elliptic flow V2 and to extract the ratio of 
the shear viscosity to the entropy density, rj/s, from sim- 
ulations of Au-I-Au collisions at RHIC energy ^/s = 200A 
GeV. Simple parton-hadron duality allows one to com- 
pare our final V2 results with the data measured at RHIC. 
Agreement between the data and theory is found with 
Glauber-type minijets initial conditions for gluons and 
with as = 0.3 — 0.6. The rj/s ratio of the gluon plasma 
created varies between 0.15 and 0.08. Standard pQCD 
interactions alone can describe the generation of large V2 
values at RHIC. The small rj/s ratios found in the simula- 



tions indicate that the gluon plasma created behaves like 
a nearly perfect fluid. This can be understood by per- 
turbative QCD without resorting to exotic explanations 
such as the AdS/CFT conjecture: Gluon bremsstrahlung 
dominates and yields rapid thermalization, and therefore 
early pressure buildup and a small shear viscosity in the 
gluon gas. 

The authors thank G. Burau and M. Gyulassy for fruit- 
ful discussions. 
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